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penetration  dynamics  are  sensitive  to  the  physical  properties  of 
the  material  next  to  the  penetrator,  it  becomes  important  to 
account  for  the  breakup  and  fracture  of  the  rock  in  the  analysis 
techniques.  Accordingly,  material  models  were  developed  that 
included  fracture  and  a post-fracture  process  which  degrades  the 
physical  properties  to  that  of  a comminuted  rock. 
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calculational  method  is  employed,  wherein  the  penetrator  is 
treated  as  a rigid  body  to  compute  the  penetration  dynamics. 

Using  the  force  loadings  on  the  penetrator  obtained  from  this 
solution,  subsequent  analyses  may  be  made  of  the  interior  response 
of  the  penetrator. 

To  test  and  demonstrate  the  code,  pre-event  calculations  of  a 
specific  rock  penetration  field  test  were  made.  The  field  experi- 
ment consisted  of  the  vertical  firing  of  a 400  pound  steel  pene- 
trator into  a massive  welded  tuff  site  at  a velocity  of  1500 
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of  the  solution  results  are  provided,  covering  the  force  loadings 
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fractured!  regions  in  the  rock,  and  the  stress  and  motions  in  the 
rock.  Comparisons  with  experimental  results  from  four  field  tests 
are  drawn!  Generally  good  agreement  was  found  for  projectile 
deceleration,  peak  stresses  developed  in  the  rock,  and  the  extent 
of  the  fracture. 
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SECTION  I 


INTRODUCTION 

1.1  BACKGROUND  AND  TASK  OBJECTIVES 

A major  current  effort  in  the  improvement  of  earth  pene- 
trator  technology  is  directed  towards  providing  the  capability 
for  penetrating  relatively  hard  target  media  (e.g.  concrete 
and  medium- strength  rock)  at  impact  velocities  in  the  1000-5000 
ft/sec  range.  A part  of  this  effort  involves  the  develop- 
ment and  use  of  a "first  principle”  finite-difference  Lagran- 
gian  computer  code,  called  WAVE-L,  for  the  analysis  of  pene- 
tration dynamics.  Previous  programs  have  examined  the  utility 
and  limitations  of  WAVE-L  for  analysis  of  projectile  impact 
and  penetration  into  soil  media. The  objective  of  the 
task  reported  herein  has  been  to  modify  the  predictive  and 
analytical  capability  afforded  by  the  WAVE-L  code  and  associ- 
ated material  models  and  to  apply  it  to  analysis  of  projectile 
penetrations  into  medium- strength  rock  media. 

To  test  and  demonstrate  the  WAVE-L  code,  pre-event  calcu- 
lations of  a rock  penetration  field  test  were  required.  The 
specific  experiment  chosen  for  this  validation  test  was  the 
planned  firing  at  1500  ft/sec  of  a 400-lb,  5-ft  long,  6.5-in. 
dia,  ogive  nose  shape,  steel  penetrator  into  a massive  welded 
tuff  site  at  the  Tonopah  Test  Range  near  Mt.  Helen,  Nevada. 

A key  requisite  for  extending  code  capabilities  to  allow 
calculations  of  earth  penetration  into  rock  media  was  the  form- 
ulation of  a material  model  that  incorporates  the  breakup  and 
fracture  characteristics  occurring  in  rock  during  impact  and 
penetration  by  the  projectile.  Examination  of  material  sur- 
rounding holes  left  in  rock  by  prior  penetration  tests  shows 
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that  there  is  a region  of  highly  comminuted  rock  surrounding 
the  penetrator  hole.  Beyond  the  comminuted  region  are  zones 
of  brecciated  and  sheared  rock,  with  the  degree  of  fracture 
diminishing  with  increasing  radius."*’^  A typical  fracture 
pattern  is  illustrated  in  Figure  1.  These  failure  patterns 
are  important  because  the  physical  properties  of  fractured 
rock  are  in  general  different  than  those  of  intact  rock. 

The  penetration  dynamics  will  be  influenced,  perhaps  strongly, 
by  these  altered  properties,  since  target  material  near  the 
advancing  projectile  nose  is  fractured,  rather  than  intact. 

It  is  of  particular  importance  that  the  comminuted  rock  has 
a lower  failure  surface  (shear  strength) , since  this  is  one 
of  the  most  sensitive  material  properties  influencing  pene- 
tration dynamics.®  Accordingly,  it  was  therefore  a necessary 
preliminary  objective  of  the  task  reported  herein  to  develop 
constitutive  relations  which  allow  assessment  of  the  extent 
of  failure,  which  has  occuri'ed  at  a given  time  and  location 
in  the  code  solution  of  penetration,  and  which  uses  this 
information  to  degrade  the  rock  properties  in  an  appropriate 
manner.  The  bases  for  the  overall  material  model  for  both  the 
intact  and  fractured  rock  were  data  from  in  situ  seismic  sur- 
veys and  data  from  a series  of  laboratory  tests  which  measured 
the  constitutive  properties  of  welded  tuff  samples  obtained 
from  the  Tonopah  site.^’® 

In  this  task,  the  specific  information  which  was  sought 
from  the  welded  tuff  penetration  analysis  was  the  deceleration 
history  and  the  loads  on  the  projectile,  and  the  media  response 
during  the  early  stage  of  the  penetration  event.  (For  nominally 
homogeneous  targets,  peak  loads  on  a penetrator  are  usually 
reached  by  the  time  the  nose  section  is  fully  embedded  in  the 
ground.)  These  results  could  then  be  compared  with  measured 
data  during  the  field  test. 
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1.2  TECHNICAL  APPROACH 


Through  experience  with  previous  earth  penetration  in- 
vestigations, it  has  been  found  to  be  efficient  to  decouple 
the  problem  into  separate  analyses  of  (1)  the  overall  pene- 
trator  dynamics  and  media  response  and  (2)  the  interior  re- 
sponse of  the  penetrator  itself.  This  approach  is  shown  in 
Figure  2.  The  key  assumption  in  this  technique  is  that  the 
penetrator  behaves  essentially  as  a rigid  body  for  the  pur- 
poses of  penetration  (Phase  I)  analyses.  (Examination  of  the 
projectile  following  the  firing  into  welded  tuff  indicated 
that  no  significant  plastic  deformation  occurred;  the  results 
of  this  study  and  other  rock  penetration  calculations  using 
the  decoupled  technique  have  shown  reasonable  agreement  with 
experiment,  indicating  that  the  rigid-body  assumption  does 
not  introduce  any  major  errors  in  the  analysis,  so  long  as 
the  deformations  of  the  projectile  are  confined  to  small, 
elastic  distortions.  If  the  impact  conditions  or  projectile 
characteristics  are  such  that  large  projectile  deformations 
occur,  deformable -body  penetration  solutions  must  be  performed 
in  place  of  using  ‘he  decoupled  technique.)  Rigid-body  pene- 
tration analyses  are  performed  with  WAVE-L,  a finite-difference, 
Lagrangian  code.  For  interior  response  (Phase  II)  analyses, 
either  finite-difference  or  finite-element  codes  may  be  used. 
Additional  discussion  of  this  decoupled  approach  may  be  found 
in  References  3 and  11. 

The  investigation  reported  herein  is  confined  to  Phase  I, 
viz. , rigid-body  penetration  into  rock  media.  The  output  from 
the  Phase  I calculations  presented  herein  could  be  used  as 
boundary  conditions  in  an  analysis  of  the  projectile’s  interior 
response. 

In  the  penetration  field  test,  which  employs  a vertical 
gun,  the  projected  mass  includes  a sabot  which  remains  attached 


PHASE  I: 


Analysis  of  penetration  dynamics 


Penetrator  Is  ass'imed  to  be  a rigid  body 
having  same  mass  and  shape  as  the  actual 
projectile. 

This  assumption  Is  justified  when  the 
distortions  of  the  projectile  during 
penetration  are  elastic  and  so  small 
that  they  do  not  significantly  affect 
the  penetration  dynamics. 


From  the  Phase 
the  projectile 
are  determined 
position  along 


I analysis,  the  forces  on 
surface  during  penetration 
(as  functions  of  time  and 
the  surface). 


PHASE  II: 


V f3(t) 


fjCt) 


Analysis  of  dynamic  structural  response 


Penetrator  Is  now  assumed  to  be  a deform- 
able body.  In  accordance  with  the  actual 
configuration  and  properties  of  Its 
materials.  The  response  of  the  projectile 
and  of  Its  internal  components,  as  they 
are  subjected  to  the  surface  force 
loading  from  Phase  I,  is  determined. 


Figure  2.  Decoupled  Approach  for  Making  Analyses  of  Penetration  and 
Response  Dynamics  of  Earth  Penetrators 
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to  the  projectile  until  the  sabot  impacts  the  ground,  i.e., 
during  the  first  5 ft  of  penetration.  The  projectile  weighed 
400  lb  and  the  sabot  weighed  117  lb.  The  total  (517  lb)  was 
used  as  the  penetrator  weight  in  the  calculations.  Since  the 
calculations  covered  just  the  early  stages  of  penetration, 
the  reduction  in  weight  occurring  when  the  sabot  is  stripped 
off  did  not  need  to  be  accounted  for. 

1.3  SUMMARY  OF  RESULTS 

Work  on  this  task  was  conducted  in  four  segments,  as 
summarized  in  the  following: 

a.  A new  material  model  for  rock  media  was  formulated 
which  includes  a fracture  criterion  and  post  - fracture  degrad- 
ation of  properties.  A specific  constitutive  relation  for 
welded  tuff  was  then  determined.  This  relation  was  based 

on  in  situ  seismic  and  laboratory  data.  However,  initially 
the  only  available  data  on  the  properties  of  crushed  tuff 
were  failure  data  at  low  pressures  (P < 5000  psi).®  The 
failure  envelope  at  higher  stresses  had  to  be  estimated  by 
extrapolation.  This  led  to  an  initial  relationship  with 
relatively  weak  properties  for  crushed  tuff. 

b.  Using  the  new  material  model  developed  in  (a),  an 
exploratory  penetration  calculation  was  conducted  of  the 
full-scale  penetrator  firing  into  the  welded  tuff  Tonopah  site. 
To  provide  a shorter  run-up  time  to  peak  loadings,  this  solu- 
tion was  started  with  the  nose  tip  of  the  penetrator  buried 

5 inches  into  the  ground.  This  exploratory  solution  was  carried 
out  to  a time  of  .64  msec  and  a penetration  of  16  in.  The 
distortional  field  from  the  solution  at  this  time  is  depicted 
in  Figure  3.  Regions  of  totally  shattered,  severely  fractured, 
lightly  fractured,  and  intact  material  are  indicated  in  the 
figure.  These  regions  are  seen  to  generally  correspond  with 
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Figure  3.  Lagrangian  Grid  and  Fracture  Pattern  in  Welded  Tuff, 
Exploratory  Solution 
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the  observations  of  the  rock  around  the  holes  from  prior 
penetration  tests  in  welded  tuff  (Figure  1).  These  results 
indicate  that  the  material  model  provides  the  correct  quali- 
tative description  of  the  rock  fracture  during  penetration. 

The  exploratory  solution  also  shows  that  the  penetrator 
is  enveloped  by  crushed  tuff  material  during  the  penetration 
process.  This  confirms  the  necessity  of  incorporating  the 
degraded  properties  of  the  comminuted  material  into  the  model 
(including,  in  particular,  its  strength  and  frictional  char- 
acteristics, since  these  factors  bear  heavily  on  the  pene- 
tration dynamics).  Due  to  the  potential  importance  of  the 
post-fracture  properties  seen  in  this  preliminary  solution, 
additional  laboratory  tests  on  crushed  tuff  were  requested. 
Higher  pressure  failure  data  were  subsequently  obtained  by 
Terra  Tek*^  and  the  Waterways  Experiment  Station  (WES).*^ 

c.  These  data  show  that  the  failure  surface  for  crushed 
material  under  high  confining  pressures  is  only  15-20%  below 
the  failure  surface  of  intact  material.  A revised  constitutive 
relationship  was  therefore  formulated  to  match  the  new  experi- 
mental determination  of  the  failure  envelope  of  crushed  tuff. 
The  hydrostat  for  intact  tuff  was  also  re-fit  at  this  time 

to  reflect  a better  interpretation  of  the  experimental  data. 

In  addition,  the  friction  rule  governing  the  shearing  stresses 
applied  to  the  penetrator  surface  was  modified  to  conform  more 
directly  with  test  data  obtained  by  Systems,  Science,  and 
Software  on  the  frictional  loads  developed  by  a steel  wheel 
sliding  over  welded  tuff  samples. 

d.  The  principal  penetration  calculation  was  then  con- 
ducted, using  the  revised  constitutive  relationship  for  welded 
tuff.  The  solution  was  carried  out  to  a time  of  1.3  msec  and 
a penetration  depth  of  23  in.  By  this  time,  the  19.5-in.  nose 
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was  completely  buried,  and  the  projectile  deceleration  had 
peaked  and  was  dropping. 


I 


1.3.1  Computational  Results 

Figure  4 shows  the  time  histories  of  the  axial  forces 
on  the  penetrator  and  the  rigid-body  deceleration  experienced. 
Friction  accounted  for  35-45%  of  the  total  deceleration.  (The 
periodic  "blips"  in  the  curves  coincide  with  rezones  of  the 
computational  grid.  These  are  performed  to  correct  severe 
distortions  of  cells  near  the  penetrator  as  they  are  dragged 
downward  by  the  strong  frictional  forces.) 

Time  histories  of  the  penetrator  velocity,  depth  of 
penetration,  and  engaged  length  (the  vertical  dimension  of 
the  "wetted"  surface  of  contact  between  the  penetrator  and 
the  rock)  are  shown  in  Figure  5.  As  is  usually  observed  during 
high-velocity  penetration,  the  steady  value  of  engaged  length 
(17.9  in.)  is  less  than  the  nose  length  (19.5  in.).  Overall 
results  of  the  solution  are  summarized  in  Table  1.  At  the 
end  of  the  solution  at  1.3  msec,  the  penetrator  velocity  is 
dropping  rapidly  and  the  stress  levels  in  the  tuff  are  de- 
clining; there  is,  therefore,  no  reason  to  expect  that  the 
forces  on  the  penetrator  will  increase  at  a later  time.  (For 
high-g  penetrations,  after  the  nose  is  fully  engaged  and  the 
initial  deceleration  peak  has  been  passed,  the  decrease  in 
force  due  to  the  rapidly  declining  projectile  velocity  will 
generally  outweigh  any  increases  in  force  due  to  increased 
confinement  at  greater  penetration  depths.  A sudden  increase 
in  deceleration  will  occur,  of  course,  when  the  mass  of  the 
sabot  is  stripped  away  at  about  5 ft  depth.)  Figure  6,  which 
shows  the  Lagrangian  grid  at  1.3  msec,  delineates  the  regions 
of  comminuted,  severely  fractured,  lightly  fractured,  and 
intact  material.  Tlie  velocity  and  principal  stress  fields 
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Figure  4 


Axial  Forces  on  Penetrator  and  Deceleration  Histories 
from  Principal  Calculation 


Axial  Forces,  kips 


VELOCITY  FT/SEC 


TIME  MSEC 


Figure  5.  Time  Histories  of  Penetrator  Velocity,  Depth  of  Penetration, 
and  Engaged  Length  from  Principal  Calculation 
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Depth  of  Penetration  or  Engaged  Length,  in. 
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Figure  6.  Lagranglan  Grid  and  Fracture  Pattern 
in  Welded  Tuff,  Principal  Solution 
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TABLE  1.  SUMMARY  OF  RESULTS  FROM  PRINCIPAL  CODE  SOLUTION 
OF  PENETRATION  INTO  WELDED  TUFF 


Peak  penetrator  deceleration 

Peak  normal  stress  on  penetrator 

Peak  tangential  stress  on 
penetrator 

Peak  axial  force/unit  length 

Peak  radial  force/unit  angle/ 
unit  length 

Frictional  force  contribution 
to  deceleration 

Engagement  (wetted)  length 
(nose  tip  to  separation) 

Peak  compressive  stress  in  tuff 

Range  of  comminution 

Range  of  cracking  on  surface 
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5200  g's  at  1.15  msec 
6.8  kb  at  1.08  msec 
1.04  kb  at  1.08  msec 

260  kips/in. , 

5 in.  back  from  tip 

110  kips/rad/in. , 

7.5  in.  back  from  tip 

35-45% 

17.9  in. 

10  kb 

1 ft  radius  (at  1.3  msec) 
3.3  ft  radius  (at  1.3  msec) 


] I 


I 

i 

I , 

■■  1 
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at  this  time  are  shown  in  Figures  7 and  8.  The  "dragdown" 
of  the  tuff  near  the  penetrator  surface,  caused  by  strong 
friction,  is  evident  in  Figures  6 and  7. 


1.3.2  Experimental  Comparisons 

Subsequent  to  the  above  code  calculations,  the  four  test 
firings  listed  below  were  conducted  by  Sandia  Corporation  at 
the  Tonopah  site:*® 

Penetration  Field  Tests  in  TTR  Welded  Tuff 


Shot  No. 

Test  Date 

Proj  ectile 

Penetration 

Velocity 

Depth 

1 

1 Jul  1975 

1644  ft/sec 

11.0  ft 

2 

19  Aug  1975 

1640  ft/sec 

11.0  ft 

3 

25  Feb  1976 

1567  ft/sec 

11.57  ft 

4 

28  Apr  1976 

1350  ft/sec 

8.6  ft 

No  telemetry  data  on  projectile  deceleration  were 
received  from  the  first  three  tests  due  to  various  electronic 
malfunctions.  Projectile  deceleration  records  from  both  a 
piezo- resistive  gage  and  a quartz  gage  were  obtained  from 
the  last,  1350  ft/sec,  test.*®  These  data,  in  terms  of  decel- 
eration vs  depth  of  penetration,  are  shown  in  Figure  9.  The 
piezo-resistive  gage  apparently  was  ringing;  the  quartz  gage 
produced  a comparatively  smooth  record  and  is  considered  by  Sandia 
experimenters  to  be  the  more  representative.  The  first  peak 
decelerations  are  3680  g's  for  the  quartz  gage  and  4370  g's  for 
the  piezo- resistive  gage.  The  subsequent  peak  at  5 ft  depth  is 
presumably  associated  with  the  stripping  off  of  the  sabot.  For 
comparison,  the  code  solution  results,  for  an  impact  velocity 
of  1500  ft/sec,  are  also  shown  on  the  figure.  A rough  scale- 
down  of  the  code  results  to  the  experimental  impact  velocity 
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500  fc/sec 


Figure  7.  Particle  Velocity  Field  in  Welded  Tuff, 
Principal  Solution 
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Figure  8. 


Principal  Stress  Field  in  Welded  Tuff, 
Principal  Solution 
22 
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Figure  9.  Comparison  of  Computed  and  Experimental  Projectile 
Deceleration  vs  Penetration  Depth 
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of  1350  ft/sec  is  also  plotted,  where  the  deceleration  was 
simply  scaled  directly  with  the  impact  velocity,  giving  a 
reduction  of  10%.  This  comparison  shows  that  the  code  solu- 
tion is  giving  a reasonable  prediction  during  the  first  stage 
of  penetration,  up  to  and  beyond  the  first  deceleration  peak. 

It  would  be  of  interest  to  carry  out  the  code  solution  to 
greater  penetration  depths  to  permit  additional  comparisons. 

Stress-time  waveforms  were  recorded  for  all  four  firings 
from  gages  placed  in  the  rock  at  depths  of  12  and  36  in. , and 
at  horizontal  ranges  from  10  to  65  in.  The  measured  peak 
amplitudes  of  the  shock  from  the  penetrating  projectile  vs 
horizontal  range  are  shown  in  Figure  10.  For  comparison,  the 
code  solution  results,  for  an  impact  velocity  of  1500  ft/sec, 
of  peak  compressive  stress  vs  range  at  a depth  of  12  in.  are 
also  shown  in  the  figure.  The  experimental  data  bracket  the 
code  results,  but  the  data  band  is  so  broad  that  the  comparison 
is  not  very  significant. 

A comparison  of  experimental  and  calculated  stress-time 
histories  at  a depth  of  12  in.  is  shown  in  Figure  11.  Again, 
there  is  such  a large  variability  in  the  experimental  waveforms 
that  it  is  difficult  to  evaluate  the  code  results  on  this  basis 
There  are  three  signals  that  have  been  tentatively  identified 
by  Sandia  in  some  of  the  experimental  records:  (1)  a low 
amplitude  disturbance  from  the  initial  impact, (2)  the  main  wave 
called  the  "crushing”  signal,  from  the  penetrating  projectile, 
and  (3)  a later-time  wave  emanating  from  the  impact  of  the 
sabot.  As  seen  in  Figure  11,  these  separate  waves  are  not 
readily  distinguishable  in  most  of  the  waveforms;  some  of  the 
"signals"  may  be  attributable  to  cable  noise.  Three  mean- 
stress  time  histories  are  shown  from  the  code  solution,  at 
horizontal  ranges  of  12,  18,  and  24  in.  These  show  both  the 


24 


Computed,  12  depth.  Impact  Velocity  = 1500  ft/sec 
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Experimental  Data  Deslenatlon 


Shot •I-IA 


A - 12"  depth 
B " 36"  depth 


« Shots  1,2,3,  J 

Impact  ' 

Velocity  “ 1567-16A4  ft/sec 


* Shot  4,  Impact 

Velocity  • 1350  ft/scc 
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Figure  10.  Comparison  of  Computed  and  Experimental  Peak 
Stress  vs.  Horizontal  Range 


precursor  and  main  waves.  The  main  wave  occurs  at  an  earlier 
time  than  that  shown  in  most  of  the  experimental  data. 

A post-test  survey  of  the  crater  region  from  the  first 
penetration  test  has  been  reported  by  WES.*^  The  diameter  of 
the  crater  at  the  surface  is  5.8  ft,  with  surface  cracks  ex- 
tending a few  inches  beyond.  This  is  in  good  agreement  with 
the  6. 6- ft  diameter  of  the  region  of  cracked  material  on  the 
surface  obtained  from  the  calculation.  In  the  experiment,  the 
tuff  next  to  the  penetrator  is  highly  comminuted,  ranging  from 
"powder  to  fractional -inch  fragments,  and  is  very  tightly 
packed."*^  Samples  of  comminuted  material  were  taken  at 
ranges  up  to  14  in.  from  the  penetrator.  Samples  taken  7 in. 
or  closer  to  the  penetrator  indicated  that  60-95%  of  the  mat- 
erial consists  of  sand-size  (<5  mm)  particles.  A sample  1 ft 
from  the  penetrator  contained  about  30%  sand-size  particles. 
The  calculation  (Figure  6)  predicted  that  highly  comminuted 
material  would  extend  to  about  1 ft  radius  (9  in.  behond  the 
penetrator) , with  highly  fractured  material  extending  an 
additional  6 in. 

The  projectile  was  recovered  intact,  with  no  apparent 
structural  damage.  The  paint  had  been  completely  stripped  off 
and  the  nose  was  slightly  abraded.  The  projectile  surface  was 
partially  covered  with  comminuted  tuff  material.  Along  the 
hole  wall  next  to  the  penetrator  there  were  patches  of  metal 
film  and  comminuted  rock  fused  together.  This  indicates  that 
sliding  friction  at  the  steel/rock  interface  is  significant 
to  the  penetration  dynamics. 

An  examination  of  the  amount  of  wear  on  the  projectile 
nose  has  been  conducted  by  Sandia;  the  wear  profile  was  found 
to  be  in  close  agreement  with  the  shape  of  the  applied  stress 
distribution  calculated  in  the  code  solution. 
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1.3.3  Differences  Between  Conditions  in  Calculation  and 
Experiment 

In  view  of  the  scattered  and  somewhat  inconsistent  nature 
of  the  test  data,  direct  comparisons  between  the  calculation 
and  the  experiment  are  difficult  to  evaluate.  There  were  also 
some  unintentional  differences  between  conditions  in  the  ex- 
periments and  in  the  code  simulation  which  affect  direct  com- 
parisons. (It  was  a program  requirement  that  the  code  solution 
be  performed  prior  to  the  test.)  These  differences  could  be 
alleviated  by  repeating  the  solution  using  impact  conditions 
identical  to  these  experienced  in  the  test,  but  such  a repe- 
tition does  not  seem  justified  due  to  questions  about  the  data. 

a.  Tb“  code  solution  was  run  using  the  anticipated 
impact  velocity  of  1500  ft/sec.  However,  in  the  test  for 
which  deceleration  records  were  obtained  (Shot  No.  4),  the 
measured  impact  velocity  was  10%  less,  or  1350  ft/sec.  An 
adjustment  to  the  calculationa 1 results  to  obtain  an  esti- 
mate of  projectile  deceleration  for  the  lower  velocity  has 
been  made  in  Figure  9. 

b.  The  projectile  in  the  code  simulation  has  a pointed 
nose  tip,  i.e.  the  ogival  shape  of  the  nose  extended  to  the 
tip.  The  test  projectile  had  a beveled  nose  tip,  i.e.  a 
short,  relatively  blunt,  conical  shape  at  the  tip  of  the 
ogival  nose.  However,  an  ogival  pointed  cap  was  glued  on 
over  the  bevel  prior  to  the  firing.  This  cap  piece  was  not 
recovered  and  probably  broke  off  early  in  the  event.  It 
would  be  preferable  to  use  an  integral  nose  tip  in  future 
comparisons . 

c.  The  experimental  impact  angle  was  slightly  oblique 
(2®)  and  the  projectile  trajectory  in  the  ground  veered  away 
from  the  vertical  on  some  shots.  The  final  inclination  angle 


of  the  recovered  projectiles  was  as  much  as  10”.  The  veering 
of  the  projectile  path  may  have  been  due  to  the  initial  yaw 
or  to  inhomogenieties  in  the  ground.  The  initial  small  impact 
angle  probably  did  not  cause  serious  discrepancies  in  the 
measured  axial  acceleration  or  in  the  magnitudes  of  the 
stresses  measured  in  the  ground.  However,  the  veering  does 
introduce  an  uncertainty  in  the  horizontal  range  of  the  gages 
in  the  ground  from  the  projectile  trajectory.  This  uncer- 
tainty is  indicated  by  the  horizontal  bars  on  some  of  the 
data  in  Figure  10. 

d.  The  material  model  includes  a provision  for  degra- 
dation of  the  properties  of  fractured  rock,  based  on  labora- 
tory tests  of  mechanically  crushed  samples.  The  properties 
of  material  undergoing  dynamic  high-pressure  comminution 
during  the  penetration  process  may  differ.  For  example,  the 
density  of  the  mechanically  crushed  samples  was  1.36  gm/cm^, 
whereas  the  calculations  indicate  that  the  material  next  to 
the  penetrator  remains  close  to  normal  density  (1.95  gm/cm^) 
during  the  comminution  process  and  while  it  is  acting  on  the 
penetrator.  Property  tests  on  samples  of  the  dynamically 
comminuted  material  would  be  desirable. 
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SECTION  II 
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PROBLEM  CONDITIONS  AND  MODELS  I 

» 

2.1  PENETRATOR  AND  IMPACT  CONDITIONS  I 

I 

The  problem  conditions  for  the  calculation  of  the  pene-  i 

tration  test  into  welded  tuff  at  the  Tonopah  Test  Range  are 
shown  in  Figure  12.  The  penetrator  was  steel  with  a cylindrical 
aft  body  and  an  ogival  nose  section.  The  total  vehicle  weight 
was  S17  lb,  consisting  of  the  400  lb  penetrator  and  a 117  lb 
sabot.  The  sabot  remains  attached  to  the  penetrator  until  it 
impacts  the  ground,  corresponding  to  a depth  of  penetration  of 
5 ft.  The  total  517  lb  weight  was  employed  throughout  the  num- 
erical solutions  since  they  were  terminated  at  a penetration 
depth  of  -2  ft.  For  the  analyses,  the  penetrator  impacted  the 
ground  at  normal  obliquity  with  a velocity  of  1500  ft/sec  and 
was  assumed  to  be  a rigid  body. 

2.2  TARGET  MEDIA  CHARACTERISTICS  AND  MODELING 

The  target  media  was  a inassive  site  of  high  quality  welded 
tuff,  known  as  Thirsty  Canyon  welded  tuff,  located  near  Mount 
Helen  on  the  Tonopah  Test  Range  in  Nevada.  The  properties  of 
the  site  material  were  characterized  through  a joint  program 
conducted  by  WES’’®  and  Terra  Tek.’  Our  analysis  of  these  ex- 
perimental findings  led  to  the  set  of  nominal  properties  listed 
in  Table  2.  These  are  used  to  set  the  initial  (undisturbed) 
properties  in  the  code  solution  of  the  problem.  The  values  of 
these  initial  mechanical  properties  are  somewhat  tenuous,  since 
there  is  a wide  range  between  the  field  seismic  data,  ultrasonic 
laboratory  data,  and  static  and  dynamic  laboratory  test  results. 

For  the  model,  intermediate  values  were  chosen.  The  results  of' 
a related  study,®  however,  have  shown  that  the  penetration  dyn- 
amics are  faily  insensitive  to  the  value  of  the  bulk  modulus, 
so  it  is  believed  that  the  effect  of  the  uncertainty  in  initial 
mechanical  properties  is  not  serious.  Other  geological  features 
of  the  site  material  are  discussed  in  the  WES  and  Terra  Tek  test 
reports. ’ 
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TARGET  MEDIUM 
Thirsty  Canyon  Welded  Tuff 


Figure  12.  Conditions  for  Pre-event  Calculations  of 
Penetration  Test  into  TTR  Welded  Tuff 
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TABLE  2. 


NOMINAL  PROPERTIES  OF  THIRSTY  CANYON  WELDED  TUFF 


Density, 

Water  content 
Air  voids 

Dilatational  velocity,  (Cp)^ 
Shear  wave  velocity,  (Cg)^ 
Unconfined  compressive  strength. 
Bulk  modulus, 

Shear  modulus. 

Constrained  modulus, 

Young's  modulus,  E^ 

Poisson's  modulus, 


= 1.95  gm/cm^ 

= 1.5%  (by  dry  weight) 

= 22%  (by  volume) 

= 2864  m/sec  (9398  £t/sec) 
= 1754  m/sec  (5755  ft/sec) 
Yq  = .63  kb 
= 80  kb 
= 60  kb 
= 160  kb 
= 144  kb 
= .2 


A material  model  for  use  in  code  calculations  was  developed 
for  the  welded  tuff  based  on  the  laboratory  constitutive  prop- 
erties data.  A hysteretic,  elastic-plastic,  strain-hardening 
model  with  a non-associated  flow  rule  was  used.  The  mechanical 
properties  in  loading  and  unloading  depend  on  the  current  volu- 
metric strain  and  maximum  volumetric  strain  which  the  material 
has  experienced. 

The  failure/fracture  model  used  to  simulate  the  observed 
non-linear  response  is  depicted  in  Figure  13.  A basic  yield 
surface,  curve(l),  was  set  along  the  line  in  stress  space  where 
the  onset  of  non-linear  response  in  intact  welded  tuff  was  ob- 
served in  triaxial  compression  tests.  The  increase  in  strength 
with  increasing  loading  was  modeled  by  means  of  a strain-hard- 
ening model,  which  raises  the  yield  surface  (as  a function  of 
generalized  plastic  strain)  until  a fracture  surface,  curve  dX 
or  maximum  yield  surface,  corresponding  to  the  line  in  stress  space 
where  sample  failure  occurs , is  reached.  Material  reaching  this  sur- 
face is  assumed  to  fracture.  A post-fracture  model  is  then  used 
that  gradually  degrades  the  yield  surface  as  a function  of 


exaess  generalized  plastic  strain;  i.e.,  that  amount  of  general- 
ized plastic  strain  occurring  after  the  onset  of  fracture.  The  post- 
fracture yield  surface  decreases  until  it  reaches  the  minimum 
yield  surface,  curve  (3!X  used  to  represent  crushed  welded  tuff 
in  the  principal  calculation.  (The  relatively  weak  minimum  yield 
surface,  curve  (|X  used  in  the  exploratory  penetration  calculation, 
is  also  indicated  in  the  figure.) 

The  load-unload  response  of  welded  tuff  under  isotropic 
compression  is  shown  in  Figure  14.  Load-unload  stress  paths 
for  uniaxial  strain  conditions  are  shown  in  Figure  15.  Both 
the  experimental  data  and  the  response  predicted  by  the  model 
are  shown  in  the  figures. 
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Figure  13.  Yield  Surfaces  for  Welded  Tuff  Models 


Figure  15.  Coirpartson  of  Model  with  Laboratory  Uniaxial  Strain 
Stress  Paths  for  Welded  Tuff 
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A detailed  description  of  the  material  model  is  given 
in  Appendix  A. 

A friction  rule  was  adopted  that  was  partially  based  on 

an  extrapolation  of  experimental  data  obtained  at  Systems, 

Science,  and  Software  (S^)  for  the  sliding  friction  between 

steel  and  welded  tuff  samples.^**  The  S^  tests  covered  only 

a limited  range  of  low  normal  stresses,  from  .14  to  .55  kb, 

and  were  at  sliding  velocities  less  than  100  ft/sec. 

Since  friction  cannot  exceed  the  shear  strength  of  the  rock, 

the  applied  frictional  stress  at  any  point  was  also  limited 

[ by  the  octahedral  shear  stress  associated  with  the  concurrent 

i value  of  the  yield  surface,  W = 1^2  at  that  point. 

[ v3  max  » 

I The  frictional  stress  rules  used  were: 

I 

I 

I T = Min  (Tq  + 

where  = 30  bars 

, P£  = .15 

[ These  rules  are  plotted  in  Figure  16.  Since  nearly  all 

1 the  rock  next  to  the  penetrator  is  in  the  crushed  state,  the 

i shear  stress  associated  with  the  minimum  yield  surface  is 

shown  in  the  plot.  As  may  be  seen,  the  first  friction  rule 
predominates  except  at  low  applied  stresses.  For  comparison, 
the  predominant  friction  rule  used  in  the  exploratory  pene- 
. tration  solutions  is  also  indicated. 

Based  on  the  evidence  available  now,  the  above  friction 
rule  is  believed  to  have  overspecified  the  frictional  stress 
acting  at  the  interface  in  the  calculation. 

The  calculational  results  of  the  ground  motions  near 
the  surface  (Figure  7)  appear  unrealistic  in  that  they  show  no 
debris  material  moving  upward.  Instead,  there  is  strong 
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Figure  16.  Friction  Rules  for  Welded  Tuff 


I 
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dragdown  or  sticking  of  material  in  contact  with  the  pene- 
trator.  The  results  of  the  exploratory  solutions,  where  the 
frictional  stress  was  considerably  smaller,  gave  a velocity 
pattern  which  appears  to  be  physically  more  realistic  (as 
shown  in  Figure  47,  page  78). 

Metallurgical  examination  of  the  recovered  penetrator 
has  shown  that  melting  of  a few  mils  of  metal  occurred.  In 
addition,  patches  of  metal  film  and  comminuted  rock  fused 
together  were  found  along  the  hole  wall.  The  presence  of  a 
melted  surface  layer  should  reduce  the  sliding  friction. 

Use  of  a lower  friction  rule  would  also  decrease  the 
calculated  penetrator  deceleration,  bringing  it  into  closer  ^ 

agreement  with  the  experimental  results.  ' 

i 
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SECTION  III 


PRINCIPAL  COMPUTATION  OF  PENETRATION  INTO  WELDED  TUFF 

Two  numerical  solutions  of  the  penetrator  firings  at  the 
Tonopah  site  were  conducted.  An  exploratory  calculation  was 
first  performed  to  evaluate  the  material  model  (in  particular, 
the  workings  of  the  comminution  and  post-fracture  treatments), 
and  to  provide  tentative  pre-test  information  on  expected 
penetrator  stress  loadings  and  stress  levels  in  the  ground. 

This  solution  is  summarized  in  Section  1.3;  additional  results 
are  given  in  Appendix  C.  Following  revisions  to  the  material 
model,  based  on  new  experimental  data  regarding  the  failure 
properties  of  crushed  tuff,  the  principal  penetration  calcu- 
lation was  conducted.  Results  of  this  principal  solution  are 
described  here. 

3.1  WAVE-L  CODE 

The  WAVE-L  computer  program  was  employed  for  the  penetration 
calculations.  WAVE-L  is  a two-dimensional  code  which  solves 
the  equations  of  motion  for  elast ic -plast ic  bodies,  using  a 
finite  difference  Lagrangian  technique.  The  basic  mathematical 
formulation  has  been  described  by  Wilkins. A sliding  inter- 
face, in  conjunction  with  the  friction  model,  was  used  between 
the  rigid-body  projectile  and  the  rock  media.  This  code  has 
been  previously  used  for  a number  of  rigid-body  penetration 
problems. 

3.2  COMPUTATIONAL  GRID  DESIGN 

The  initial  grid  for  the  penetration  solution  is  shown  in 
Figure  17.  The  contour  of  the  penetrator  periphery  was  resolved 
with  a total  of  50  lattice  points,  with  a spacing  in  the  vertical 
(axial)  direction  along  the  ogival  nose  of  -.4  in.  (one-eighth  of 
the  projectile  radius).  The  initial  grid  in  the  rock  media 


40 


consisted  of  ~.8  x .8  in.  cells  (one-fourth  of  the  projectile 
radius)  in  the  region  from  the  axis  of  symmetry  out  to  a radius 
of  6.5  in.  (one  projectile  diameter).  (The  results  of  a re- 
lated study  on  the  effects  of  zone  size  on  computational  results 
indicate  that  this  resolution  (zone  fineness)  sould  be  adequate.^ 
Beyond  this,  the  radial  zone  width  was  increased  in  2%  steps. 

The  vertical  zone  widths  similarly  increase  in  21  steps  below 
a depth  of  25  in.  The  grid  extended  to  a radius  of  5 ft  and 
a depth  of  6 ft.  The  entire  grid  comprised  about  3500  cells. 

As  shown  in  Figure  17,  the  grid  lines  in  the  first  two 

columns  next  to  the  axis  are  canted  upward.  As  the  projectile 

encounters  cells  and  these  are  forced  around  the  nose,  they 

tend  to  distort  downward,  particularly  when  strong  frictional 

forces  are  acting.  When  this  distortion  becomes  severe,  a local 

rezone  of  the  grid  is  made.  The  upward  canting  of  the  cells 

extends  the  time  interval  between  rezones  and  also  reduces 

1 > 2 

numerical  oscillations. 

3.3  SOLUTION  RESULTS 

The  principal  calculation  of  the  penetration  was  carried 
out  to  a time  of  1.32  msec  and  a penetration  depth  of  23.1  in., 
which  was  beyond  the  point  of  peak  loading  on  the  penetrator. 

A summary  of  key  results  from  the  solution  is  given  in  Table  1 
(page  18  in  Section  1.3). 

3.3.1  Penetrator  Deceleration  and  Loading 

The  computed  deceleration,  velocity,  and  displacement 
(depth)  histories  of  the  penetrator  appear  in  Figures  4 and  5 
(Section  1.3).  A peak  value  of  5200  g's  deceleration  was  reached 
at  1.15  msec,  at  a penetration  depth  of  20.3  in.  The  axial 
forces  applied  to  the  penetrator  are  also  shown  in  Figure  4. 
Frictional  force  along  the  penetrator/rock  interface  provided 
an  increasing  proportion  of  the  total  decelerating  force,  from 
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33%  to  45%,  as  shown  in  Figure  18.  The  frictional  contribution 
can  be  expected  to  increase  during  the  initial  entry,  as  the 

engaged  length  increases.  A plot  of  the  engaged  length  vs  time 

* 

is  included  in  Figure  5.  The  engaged  length  reaches  a steady 

level  (~17.9  in.)  as  the  hole  wall  separates  from  the  pene- 

trator  nose.  The  percent  of  the  frictional  contribution  in  the 

total  decelerating  force  continues  to  increase  after  separation  : 

occurs,  since  the  normal  forces  near  the  penetrator  tip  are 

declining  as  the  projectile  velocity  drops.  ? 

Time  histories  of  the  radial  forces  applied  to  the  pene- 
trator are  shown  in  Figure  19.  Tlie  peak  radial  force  was  1150 
kips/radian,  occurring  at  1.15  msec. 

Distributions  of  the  normal  and  tangential  (frictional)  • 

stresses  applied  to  the  penetrator  at  several  times  during  the  \ 

solution  are  shown  in  Figure  20.  These  stresses  are  always  \ 

highest  near  the  nose  tip.  The  peak  normal  stress  was  6.8  kb. 

Associated  force  distributions  are  shown  in  Figures  21  and  22.  • 

After  the  penetrator  is  fully  engaged,  the  peak  axial  force  ] 

per  unit  length  of  the  penetrator  is  sustained  at  about  5 in. 
behind  the  penetrator  tip.  The  peak  axial  force  was  260  kips/in., 
and  the  peak  radial  force  was  110  kips/radian/ in . 

Representative  time  histories  of  the  normal  and  shear 
stresses  applied  at  nine  individual  points  on  the  penetrator 
are  shown  in  Figures  23  and  24.  The  high  frequency  oscillations  J 

in  the  stresses  applied  to  points  near  the  nose  are  numerical.  j 

They  are  due  to  the  build-up  of  stresses  in  each  new  computational  j 

cell  encountered  as  the  nose  tip  penetrates  the  grid. 

Profiles  of  the  radial  stress  in  the  rock  vs  radial  dis- 
tance (range)  from  the  penetrator  at  a time  of  1*32  msec  are 
shown  in  Figure  25.  The  dots  on  the  curves  correspond  to  the 
radial  distance  to  the  cell  centers. 
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Figure  19.  Radial  Forces  on  Penetrator  vs  Time 
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Figure  20.  Distributions  of  Normal  and  Tangential  Stress 
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Figure  2T.  Norniiil  Stress  Loadings  at  Several  Points  on  the  Penetrator  Surface 
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3.3.2  Ground  Response 


3. 3. 2.1  Peak  Stress 

The  peak  compressive  stresses  attained  in  the  welded  tuff 
are  shown  in  Figures  26  to  28.  Peak  stress  vs  depth  along 
the  central  axis  and  at  five  ranges  is  plotted  in  Figure  26. 

The  highest  stress  experienced  in  the  ground  was  9.8  kb, 
occurring  near  the  axis  at  a depth  of  14  in. , at  a time  of 
1 msec.  Peak  stress  vs  radius  at  depths  of  6 and  12  in.  is 
plotted  in  Figure  27.  A contour  plot  of  the  peak  compressive 
stress  attained  in  the  ground  as  of  1.32  msec  is  shown  in 
Figure  28. 

A contour  plot  of  the  peak  tensile  stress  experienced 
in  the  ground  by  1.32  msec  is  shown  in  Figure  29.  The  peak 
tensions  computed  probably  exceed  somewhat  the  tensile  strength 
of  welded  tuff  (measured  at  33  bars  in  simple  tension). 

3. 3. 2. 2 Stress  Profiles  and  Stress,  Velocity,  and  Density 
Contours  in  the  Ground 

Profiles  of  the  stress  components  along  the  first  column 
of  cells  (next  to  the  central  axis  or  the  penetrator  surface)  for 
several  times  are  shown  in  Figures  30  to  34. 

The  ground  response  at  1.32  msec  (at  the  end  of  the  compu- 
tation) is  illustrated  in  Figures  35  to  40,  which  are  individual 
contour  plots  of  the  stress  and  velocity  components. 

A contour  plot  of  the  density  is  shown  in  Figure  41.  The 
initial  density  of  the  welded  tuff  is  = 1.95  gm/cm^.  Only 
the  ground  material  near  the  surface  is  seen  to  have  expanded 
out  by  this  time  (1.32  msec)  in  the  event. 
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Figure  26.  Peak  Compressive  Stress  vs  Depth  for  Several  Radii 
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Figure  27.  Peak  Compressive  Stress 
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Figure  29.  Contours  of  Peak  Tensile  Stress  (bars)  Experienced 
by  1.32  msec  In  Welded  Tuff 
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Figure  31.  Stress  Profiles  along  Central  Axis  or  Penetrator  Surface, 
t • .549  msec 
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Figure  32.  Stress  Profiles  along  Central  Axis  or  Penetrator  Surface 
t * .786  msec 
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Figure  33.  Stress  Profiles  along  Central  Axis  or  Penetrator  Surface, 
t » 1.135  msec 
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Figure  34.  Stress  Profiles  along  Central  Axis  or  Penetrator  Surface 
t « 1.322  msec 
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Figure  38.  Contours  of  Shear  Stress  (kb)  in  Welded  Tuff  at  1.32  msec 
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Figure  41.  Contours  of  Density  (gm/ctn^)  In  Welded  Tuff 
at  1.32  msec 
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3. 3. 2, 3 Compression  - Extension  States 

The  property  tests  on  welded  tuff  samples  showed  that 
the  failure  surface  was  not  symmetrical  for  triaxial  compres- 
sion and  triaxial  extension  tests,  with  the  extension  failure 
surface  having  a lower  amplitude.  This  is  not  unusual  for  geo- 
logical materials;  nonetheless,  most  material  models  have  used 
symmetrical  failure  surfaces,  usually  consisting  of  a surface 
of  revolution  about  the  (1,1,1)  line  in  (o^,  02*  o^)  principal 
stress  space.  The  symmetrical  failure  surfaces  have  usually 
been  used  since  in  many  cases  only  triaxial  compression  data 
were  available , or  the  asymmetry  was  not  considered  to  be  sig- 
nificant. Considerably  more  effort  is  entailed  in  developing 
asymmetric  failure  models,  as  well,  and  data  under  torsion  or 
biaxial  loading,  which  are  required  to  determine  the  character- 
istics of  the  failure  surface  between  the  triaxial  compression 
and  triaxial  extension  lines,  are  usually  not  available. 

A general  model  has  been  worked  out  which  can  be  employed 
to  fit  data  obtained  under  arbitrary  loading  conditions.  This 
model  was  used  to  fit  Solenhofen  limestone,  for  which  a large 
body  of  data  under  conditions  of  triaxial  compression,  triaxial 
extension,  torsion,  and  biaxial  loading  was  available.  The  model 
is  described  in  Reference  20. 

For  the  welded  tuff  model,  there  was  not  time  to  invest- 
igate the  possiblity  of  using  a general  model.  However,  it 
was  decided  to  examine  the  stress  states  occurring  in  the  solu- 
tion to  see  to  what  extent  they  departed  from  the  line  of  tri- 
axial compression,  upon  which  the  symmetrical  model  was  based. 

By  using  a rotated  principal  stress  space,  where 

X = (Oj^- 202+03) /\/6 

y = C01-O3)/  'll 

z = (Oj+02  +03)/  \/3 
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and  projecting  the  resulting  states  onto  the  x,y  plane,  the 
stress  conditions  can  be  visualized,  as  shown  in  Figure  42 
(this  technique  is  more  fully  described  in  Reference  20).  If 
the  stresses  are  ordered  according  to  Oj>  02>  o^,  all  states 
will  lie  in  the  sector  shown.  The  sector  boundaries  corres- 
pond to  triaxial  extension  and  triaxial  compression  as  indi- 
cated. 

Figures  43  to  45  show  the  stress  states  occurring  at 
1.32  msec.  The  stress  states  in  each  vertical  column  of  cells 
starting  from  the  central  axis  are  plotted  (corresponding  to 
.4  to  12.2  in.  radii).  Generally,  the  stress  states  at  ground 
depths  below  the  depth  of  the  penetrator  tip  are  seen  to  lie 
close  to  the  triaxial  compression  line.  States  lying  above 
the  penetrator  tip  are  seen  to  diverge  toward  the  y axis  (line 
of  torsion).  At  the  higher  radii,  some  low-stress  extension 
states  are  seen.  Since  the  actual  yield  surface  would  have 
downward- slanting  lines  across  the  sector  from  right  to  left 
(see  Figure  42) , the  stress  states  lying  off  the  triaxial  com- 
pression boundary  would  be  forced  downward  from  their  present 
position,  limiting  the  buildup  of  deviatoric  stress.  The  mat- 
erial strengtft  is  thus  overestimated  with  the  symmetrical  model. 

3. 3. 2. 4 Energy  Partition 

The  energy  partitioning  between  the  penetrator  and  the 
ground  during  the  solution  is  given  in  Figure  46. 

Additional  computational  results  are  contained  in  Appendix 
B;  included  are  time  histories  of  stress  and  velocity  and  stress 
paths  at  several  locations  in  the  ground,  and  additional  defor- 
mation, velocity,  and  stress  field  plots  at  intermediate  times. 
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Figure  45.  Stress  States  in  Rotated  Principal  Stress  Space,  Columns  11  to  15 


3.4  Sensitivity  of  Penetration  to  Post-Fracture  Properties 

A comparison  of  the  exploratory  and  principal  solutions 
provides  an  assessment  of  the  effect  of  differences  in  the 
post- fracture  properties  of  the  tuff  on  the  penetrator  dynamics. 
As  pointed  out  previously,  the  fractured  rock  was  assumed  to  be 
considerably  weaker  for  the  exploratory  solution.  Since  all 
the  rock  next  to  the  penetrator  surface  is  in  a highly  fractured 
state,  the  properties  of  the  comminuted  rock,  rather  than  these 
of  the  intact  rock,  predominate  in  the  determination  of  the  local 
forces  on  the  penetrator.  The  frictional  stress  must  also  be 
that  appropriate  for  steel  sliding  on  comminuted  rock.  At  a 
representative  stress  level  of  2 kb,  the  yield  strength  is 
.95  kb  and  the  frictional  stress  is  .08  kb  using  the  weak  model; 
this  compares  with  a yield  strength  of  2.4  kb  and  a frictional 
stress  of  .33  kb  using  the  strong  model.  This  difference  in 
strength  and  friction  resulted  in  significantly  different  pene- 
tration results,  as  compared  in  Table  3.  The  table  lists  results 
from  the  exploratory  and  principal  solutions  when  the  penetrator 
is  at  a depth  of  16  in.,  near  the  point  of  peak  projectile 
deceleration.  This  strong  sensitivity  of  penetration  mechanics 
to  the  post- fracture  target  properties  indicates  that  this  is 
an  area  that  should  receive  careful  attention  in  future  inves- 
tigations . 
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SECTION  IV 


CONCLUSIONS 

An  objective  of  this  effort  was  to  develop  and  verify 
"first-principle”  hydrocode  techniques  for  obtaining  detailed 
descriptions  of  penetration  dynamics  occurring  for  the  normal 
impact  and  penetration  of  projectiles  into  rock  media.  The 
results  described  in  this  report  indicate  that  development  of 
the  techniques  has  been  essentially  completed. 

An  important  result  of  the  effort  was  recognition  of  the 
importance  of  the  fracture  and  comminution  of  rock  occurring 
during  projectile  penetration.  Since  the  penetrator  is  envel- 
oped by  comminuted  material  during  the  penetration  process, 
it  becomes  essential  to  incorporate  a model  of  comminution  that 
predicts  fracture  and  then  ascribes  appropriate  post -fracture 
properties  to  the  broken-up  material.  This  important  part  of 
the  model  appears  to  be  working  properly,  with  predicted 
fracture  patterns  appearing  to  be  in  general  agreement  with 
that  observed. 

Probably  the  most  uncertain  aspect  of  the  penetration 
process  concerns  the  amount  of  frictional  force  being  applied 
to  the  penetrator.  The  severe  dragdown  of  material  next  to 
the  penetrator  and  lack  of  debris  blowoff  near  the  surface 
(as  shown  in  Figure  7)  indicates  that  friction  was  probably 
overspecified  in  this  solution.  The  results  of  the  exploratory 
solution,  which  used  a model  with  a smaller  amount  of  friction, 
gave  a velocity  pattern  which  appears  to  be  physically  more 
realistic  (Figure  47).  Post-test  surveys  and  analysis  of  the 
field  event  records  should  provide  some  guidance  on  this  problem. 
It  is  known  that  considerable  heating  and  some  melting  takes 
place  at  the  projectile  surface,  since  patches  of  metal  film  and 
comminuted  rock  fused  together  were  found  along  the  hole  wall. 


OLinmin  ucmkh  and  TEotaLBOT  mwe-l  code 

RUN  NO.  2000-201.  PEMETHITION  INTO  NELOEO  TUFF  600  FT/3EC 

in  CTCIE  787  ' 


Figure  47.  Particle  Velocity  Field  in  Welded  Tuff,  Preliminary  Solution  | 
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The  accuracy  of  the  code  solutions  has  been  partially  con- 
firmed through  comparisons  with  field  test  results.  To  more 
fully  verify  the  code  solutions,  longer-durat ion  penetration 
solutions  should  be  conducted  for  actual  field  test  conditions 
for  penetration  tests  where  reliable  projectile  and  ground 
response  data  have  been  obtained. 
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APPENDIX  A 


WELDED  TUFF  MATERIAL  MODEL 

This  appendix  gives  the  mathematical  description  of  the 
material  model  developed  for  welded  tuff.  General  features 
of  the  model  are  described  in  Section  2.2.  The  model  was 
based  on  experimental  constitutive  properties  data  determined 
by  WES  and  Terra  Tek  (Refs.  A-1  to  A-3).  Differences  between 
this  version  of  the  model  and  the  model  used  for  the  explor- 
atory penetration  solution  are  indicated. 

Symbol  Definitions 

= bulk  modulus  in  loading 
By  = bulk  modulus  in  unloading,  reloading 
= shear  modulus  in  loading 
Gy  = shear  modulus  in  unloading,  reloading 
= second  invariant  of  deviatoric  stresses 
P^  = pressure  in  loading 

P = maximum  pressure  reached  by  a material  element 
Py  = pressure  in  unloading,  reloading 
Y = yield  strength,  basic  yield  surface 
degraded  value  of  yield  strength 

yield  strength,  maximum  yield  surface  (fracture  surface) 
yield  strength,  yield  surface  of  crushed  material 
strain-hardened  value  of  yield  strength 
generalized  plastic  strain 
excess  generalized  plastic  strain 

plastic  strain  component 

- 1 = natural  volumetric  strain 
maximum  volumetric  strain  reached  by  a material  element 
Poisson's  ratio  in  loading 
Poisson's  ratio  in  unloading,  reloading 
density 
normal  density 
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Hydrostat 

The  model  hydrostat  was  based  on  load-unload  paths  mea- 
sured in  isotropic  compression  tests  conducted  by  Terra  Tek 
and  WES.  This  hydrostat  differed  slightly  from  the  one  used 
in  the  first  model  due  to  a re- interpretation  of  the  data  by 
the  testing  laboratories.  The  principal  changes  from  that 
used  before  were:  (a)  the  Terra  Tek  data  was  not  shifted  to 
eliminate  an  initial  toe  and  (b)  the  representation  given  in 

I Ref.  A-1  was  used  solely  instead  of  the  previously  recommended 

set  of  hydrostatic  load-unload  stress- strain  curves  (Ref.  A-2). 

Hydrostat  Equations  - Loading 

I 

! 

! < 


Po  = aiP 


y < 0 


P„  = 


P„  = 


Pj  + a^(y-y-)  + b-Cy-y^)^  + c^(y-y^)^. 


yi<  y £ yj^  + 1,  i = 1,2, 3, 4 


The  values  of  the  constants  are: 

a = .08  Mb 
1 

b = .261408  Mb 
1 

c = 29.8702  Mb 

a^  = .0978634  Mb 
2 

b = 1.29193  Mb 
2 

c = -82.8818  Mb 
2 

a = .0989563  Mb 
3 

b = -1.18209  Mb 
3 

c = 15.0208  Mb 

3 

a = . 068  Mb 

4 

b = -.910870  Mb 
4 

c = 6.34517  Mb 
4 

K = .025  Mb 


m 

Pi 

P2 

P3 

P4 

Pc 

Pi 

P2 

P3 

P4 


.6  Mb 

0. 

.001  Mb 
.00202  Mb 
.004  Mb 
.006  Mb 

0. 

.0115 
.02145 
. 0466 

. 1 
.1 
.8 
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P5 

Pc 

y* 


Hydrostat  - Unloading,  Reloading  ) 

Pu  ■ < (“■I's) 


Pu  ■ ^ (“■I's) 

■“u  ■ ■'o  (“■Ms)  * A(u-u  )“ 


where 


A = 


a = 


_ ^max 


^^max"’^s^ 

^^max~ 

^max”*^o  ^^max' 


M < p, 


vj  > y. 


K = K 
m uo 


•'m  = 


y_  = 


u - u 

’^max 


i . y 

max  a 

Ph  - Po 
b a 


(r  y ) 
max-* 


y<s  = r y 


max 


u < u 
•^max  - *^0 


^max  ^ 


y < u 

max 


y,s  y y. 

a ^max 


y > u 
max  — b 


The  values  of  the  constants  are; 


, 1 

0 

.03 

Mb 

= 

.01 

m 

.6 

Mb 

% = 

.04 

uo“ 

.16 

Mb 

"S  = 

.8 

= 

.15 
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Hydrostatic  tension  was  limited  by  imposing  a minimum 

value  of  pressure:  = -33  bars. 

min 

Shear  Modulus  - Loading 


3B  (l-2v  ) 

Min  — — , G, 

. 2{l+v^)  ' 


where 


= (u+1) 


Shear  Modulus-Unloading,  Reloading 

luo  X 


where 


= .5  - .32  exp(-10 


Yield/Fracture  Surfaces 

The  yield  surface,  Y(P),  defines  the  limit  of  elastic 
states;  i.e.,  v/sjJ  £ Y.  For  tentative  stress  states  lying 
outside  the  yield  surface,  plastic  flow  is  computed  in  ac- 
cordance with  the  Prandtl-Reuss  (non-associated)  flow  rule. 

The  yield  surfaces  developed  for  the  welded  tuff  model 
are  shown  in  Figure  A-1. 
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Maximum  Yield  Surface  (Fracture  Surface) : 


+ ^2 


max 


= Y.  + a.(P-Pp  + bj(P-P.)"  + c-CP-P.)-^, 


Pi<P  1 Pi.i.  i = 1.  2.  3 


%ax  = ^4  ^ 


P > P; 


Basic  Yield  Surface 


Y 

Y 


= Y 


max 


P 1 Pq 


Yi  + d.(P-P.)  + e.  (P-P^)^  + £i(P-Pi)^ 


Y = Y, 


P.<P<P.  i=l  2 3 


p > p. 


The  minimum  yield  surface  was  raised  from  that  used  in 
the  first  model  to  match  new  higher-pressure  data  obtained 
for  mechanically-crushed  tuff  (Ref.  A-4  and  A-5).  The  par- 
ticle size  distribution  of  the  fragments  used  in  the  test 
samples  approximated  a coarse  sand. 


Minimum  Yield  Surface  (completely  degraded  following  fracture) 


Y . =0 

min 


-^2 


P < 0 
,3 


Vn  = - u.(P-Pi)  ^ v.(P-P.)‘-  . w.(P-P.) 


Pi<  P 1 Pi.l*  2 
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The  values  of  the  constants,  for  Y and  P in  Mb , are: 


a^  = 13. 

= -3.1376  X 10' 

= 3.13561  X 10' 

^2  ~ 2 . 53589 

b2  = -3334.43 

C2  = 4.80632  x 10^ 

a^  = 1.77131 

63  = -306.452 

C3  = -10765.7 

a.  = 1.1464 


-3.3  X 10 

0. 

2.1  X 10 
.00118 

0. 


-5 


-4 


= 2 X 10 


-4 


= 6.3  X 10 


-4 


= 6.9725  X 10 

■4 
■4 


-4 


= 

1.2 

P 

0 

= 

2.1  X 

®1 

-466.648 

Pi 

2.1  X 

fl  = 

128256. 

P2 

* 

. 001 

^2  = 

.702829 

P3 

= 

.003 

®2  = 

-162.681 

P4 

= 

.008 

^2 

18133.3 

Yo 

= 

6.3  X 

II 

. 269704 

Yi 

= 

6.3  X 

®3  ■ 

-53. 8815 

^2 

= 

.00135 

Lti 

II 

3588.15 

Y3 

.00225 

= 

.0027 

c 

II 

1.6 

p; 

0. 

^1  = 

-319.097 

1 

Wi  = 

61411.9 

P2 

. 00171 

^^2  = 

1.04741 

P3 

= 

.01 

II 

> 

-4.05443 

Y' 

_ 

0. 

«2  " 

-974.364 

0 

.-4 


Y,  = 


y:  = 


.00211 

.0067 
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...  . 


The  basic  yield  surface  strain  hardens  in  accordance  with: 


''sh  ■ ''o  * [''-'fj  [5-2  «P(-8eP]  eP>  0,  P>P„ 


where 

eP  =/(|  deP.deP.)*' 

Strain-hardened  yield  strengths  may  not  exceed  the  max- 
imum yield  surface  ’ however. 

A state  reaching  the  maximum  yield  surface  signifies 
the  onset  of  fracture.  Material  which  has  thus  fractured  is 
governed  by  the  maximum  yield  surface,  which  then  degrades  as 
a function  of  excess  generalized  plastic  strain,  » accumu- 
lated from  the  onset  of  fracture; 


^dgd  = \ax(l-ll 


The  degraded  (post-fracture)  yield  surface  has  minimum  values 

defined  by  Y„.  , a yield  surface  representative  of  crushed 
' min’  ' 

welded  tuff  (Refs.  A-3  to  A-S). 

The  hydrostatic  tension  limit  is  similarly  degraded  to 

a final  value  of  P„;„=0  for  the  crushed  tuff. 

min 

Model  Results  and  Comparisons 

Plots  of  vertical  stress  vs  vertical  strain,  mean  normal 
stress  vs  volumetric  strain,  and  stress  difference  vs  mean 
normal  stress  for  uniaxial  strain  load-unload  paths  to  pressure 
maxima  of  1,  3,  and  5 kb  computed  with  the  model  are  shown  in 
Figures  A-2  to  A-4.  The  form  of  the  model  hydrostat  up  to 
40  kb  pressure  is  shown  in  Figure  A-5. 
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Mean  Normal  Stress,  kb 
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Figure  A-5.  Model  Hydrostat  for  Welded  Tuff 
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A comparison  of  the  load-unload  hydrostat  with  the  lab- 
oratory results  is  shown  in  Figure  A-6.  The  model  is  seen 
to  provide  a close  fit  to  the  isotropic  compression  data. 

Also  shown  on  this  figure  is  the  implied  load-unload  path 
derived  from  curves  given  for  a static  uniaxial  strain  test. 

This  curve  (and  Figure  8 of  Ref.  A-1)  indicates  that  the  mat- 
erial compacts  in  shear.  Modeling  of  this  aspect  of  material 
behavior  (volume  changes  due  to  shearing)  was  not  attempted 
in  this  effort.  Comparisons  of  vertical  stress  vs  vertical 
strain  and  stress  difference  vs  mean  normal  stress  for  a uni- 
axial strain  load-unload  path  are  given  in  Figures  A- 7 and 
A-8.  Since  the  model  does  not  exhibit  shear  compaction, 
the  model  stress- strain  curve  lies  somewhat  above  the  test 
data  (Figure  A-7).  The  fit  to  the  uniaxial  strain  stress 
path  is,  however,  reasonably  close  (Figure  A-8). 

Comparisons  of  the  model  with  experimental  triaxial  shear 
stress- strain  load-unload  paths  are  shown  in  Figure  A-9.  These 
comparisons  were  made  after  the  model  development,  as  a stress- 
controlled  driver  program  was  not  available  at  that  time.  In 
the  initial  assessment  of  the  constants  for  the  strain-hardening 
portion  of  the  model,  analysis  of  the  triaxial  shear  stress- 
strain  data  indicated  a certain  rate  of  strain-hardening.  This 
rate  was  subsequently  lowered  to  provide  a better  fit  to  the 
uniaxial  strain  data.  By  using  the  faster  strain-hardening 
law  (by  changing  the  value  of  the  constant  in  the  strain- 
hardening equation  on  Page  89  from  8 to  80),  an  improved  fit 
to  the  triaxial  data  is  obtained,  as  shown  in  Figure  A-10.  (The 
model  with  the  lower  value,  giving  the  fit  shown  in  Figure  A-9, 
was  used  in  the  calculation.)  Note  that  this  model  can  produce 
a reduction  in  stress  difference  before  the  maximum  strain  dif- 
ference is  reached.  This  is  achieved  through  the  degradation 
model  following  fracture. 

Friction  Rule 

The  friction  rule  employed  for  the  welded  tuff  is  des- 
cribed in  Section  2.2. 
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Figure  A-6.  Comparison  of  Hydrostatic  Load-Unload  Stress-Strain 
Curves 
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Figure  A-9*  Comparison  of  Triaxial  Test  Principal  Stress  Difference- 
Principal  Strain  Difference  Paths 


Figure  A-10.  Comparison  of  Triaxial  Test  Principal  Stress  Difference- 
Principal  Strain  Difference  Paths  with  Increased  Rate  of 
Strain  Hardening  in  Model 
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APPENDIX  B 


ADDITIONAL  COMPUTATIONAL  RESULTS  - GROUND  RESPONSE 

B.l  STRESS  AND  VELOCITY  TIME  HISTORIES  AND  STRESS  PATHS 
IN  THE  FIELD 

Time  histories  of  stress  and  velocity  were  obtained  at 
15  locations  in  the  ground.  The  following  parameters  are 
plotted: 

a.  Radial  stress 
Vertical  stress 

Hoop  stress  Oq 

Shear  stress 

b.  Radial  particle  velocity  f 

Vertical  particle  velocity  z 

Positive  stresses  are  compressive.  Positive  r is 
radially  outward.  Positive  z is  vertically  down. 

In  addition,  the  stress  paths  ( Jj2  vs  P)  at  these  sta- 
tions were  obtained.  In  these  plots,  the  dots  indicate 
loading  and  the  pluses  indicate  unloading.  Also  shown  on 
these  plots  are  the  initial  position  of  the  maximum  yield 
surface  (Y_„^)  and  the  minimum  yield  surface  (Y_-  ) . (The 
maximum  yield  surface  degrades  toward  the  minimum  yield  surface 
as  a function  of  generalized  plastic  strain.)  These  plots 
can  provide  guidance  as  to  the  type  of  lab  tests  which  should 
be  run  to  simulate  the  loadings  experienced  in  a penetration 
event.  The  loadings  are  characterized  by  a rapid  rise  in 
deviatoric  stress  at  low  mean  stress,  which  are  steeper  than 
attained  in  a triaxial  compression  test  at  zero  confining 
pressure.  Also,  the  unload  paths  maintain  a high  deviatoric 
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stress,  often  staying  on  the  yield  surface.  (Note  that  the  pop- 
ular uniaxial  strain  load-unload  lab  tests  are  not  at  all  rep- 
resentative of  the  stress  paths  indicated  by  the  penetration 
calculations.  Since  this  behavior  had  been  discovered  in  pre- 
vious soil  penetration  calculations,  the  laboratory  tests  for 
this  project  were  set  up  to  include  the  more  representative 
triaxial  tests.) 

The  locations  of  the  field  stations  are  indicated  in 
Figure  B-1.  The  plots  obtained  at  each  of  these  stations  are 
given  in  Figures  B-2  to  B-46.  Some  of  the  "choppiness"  in  the 
plots  along  the  axis  are  due  to  repositioning  of  the  points  at 
rezones . 

B-2  FIELD  PLOTS  OF  DEFORMATION , VELOCITY , AND  PRINCIPAL  STRESS 

In  addition  to  the  field  plots  shown  in  Figures  6 to  8 
depicting  the  final  solution  results,  field  plots  were  obtained 
at  intermediate  times  during  the  solution.  Field  plots  for  pene- 
tration depths  of  5,  10,  and  15  in.  are  shown  in  Figures  B-47  to 
B-55.  For  each  depth,  plots  of  the  Lagrangian  grid  and  fracture 
pattern,  velocity  vectors,  and  principal  stresses  are  shown. 

The  plots  of  the  grid  indicate  the  development  of  fracture 
in  the  target;  the  following  symbols  are  used  to  indicate  the 
state  of  each  cell  of  material; 

Fracture  Symbol 


Intact,  material  has  not  fractured  (has  not 
reached  Y^^^) 


□ 

a 

a 


Lightly  fractured  or  cracked  material. 


ex 


Severely  fractured  material,  ^ 1% 

Comminuted  material,  fractured  and  completely 

degraded  to  Y„. 

° min 
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There  is  no  particular  significance  to  the  dividing  point 
of  II  in  excess  generalized  plastic  strain,  except  to 

quantitatively  indicate  the  degree  of  fracture.  There  is 
actually  a smooth  gradation  of  , and  thus  severity  of 
fracture,  in  the  field. 

The  velocity  vector  field  plots  show  the  direction  and 
magnitude  of  the  particle  velocity  at  each  lattice  point  in 
the  computing  grid.  (For  clarity  in  viewing  the  ground  re- 
sponse, the  velocity  vectors  of  the  projectile  are  not  shown.) 

In  the  stress  field  plots,  the  principal  components  of 
the  stress  tensor  for  each  cell  are  shown,  as  follows:  The 
magnitude  of  the  two  principal  stresses  in  the  r-z  plane  are 
plotted  in  their  corresponding  principal  directions.  The 
third  principal  stress  (in  the  azimuthal  direction)  is  plotted 
along  the  line  bisecting  the  other  two  principal  directions. 
Vectors  pointing  to  the  right  are  compressive,  to  the  left, 
tensile.  An  example  of  how  a stress  tensor  is  plotted  is 
sketched  below: 
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Figure  B-1.  Locations  of  Time-History  Stations 
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Figure  B-11.  Stress  Components  at  Station  15 
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Figure  B-12.  Velocity  Components  at  Station  15 
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Figure  B-15.  Velocity  Components  at  Station  16 
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Figure  B-17.  Stress  Components  at  Station  17 
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Figure  B-18.  Velocity  Components  at  Station  17 
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Figure  B-20.  Stress  Components  at  Station  18 
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Figure  B-21.  Velocity  Components  at  Station  18 
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Figure  B-24.  Velocity  Components  at  Station  19 
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Figure  B-29.  Stress  Components  at  Station  21 
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Figure  B-33.  Velocity  Components  at  Station  22 


138 


CNLireWM  MKMGH  NM  nOMBLan,  ilC. 

WM  MS.  no-age.  rtmmtim  into  mmo  nvr 
smTiBN  aa 


0 O.J  0.2  0.3  0.11  O.S 

PRESSURE  KB 


; Figure  B-34.  Stress  Path  vs  P)  at  Station  22 

\ 

i 


139 


0 0.5  1.0  1.5 


TIME  nSEC 

Figure  B-35.  Stress  Components  at  Station  23 
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Figure  B-36.  Velocity  Components  at  Station  23 
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Figure  B-39.  Velocity  Components  at  Station  2A 
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Figure  B-40.  Stress  Path  ***■  24 
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Figure  B-43.  Stress  Path  vs  P)  at  Station  25 
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Figure  B-A5.  Velocity  Components  at  Station  26 
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Figure  B-46.  Stess  Path 
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Figure  B-A7.  Lagranglan  Grid  and  Fracture  Pattern  In  Welded  Tuff 
at  5 In.  Depth  of  Penetration 
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Figure  B-48.  Particle  Velocity  Field  in  Welded  Tuff 
at  5 in.  Depth  of  Penetration 
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Figure  B-50.  Lagrangrlan  Grid  and  Fracture  Pattern  In 

Welded  Tuff  at  10  In.  Depth  of  Penetration 
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Figure  B-51.  Particle  Velocity  Field  in  Welded  Tuff  at 


10  in.  Depth  of  Penetration 
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Figure  B-52.  Principal  Stress  Field  in  Welded  Tuff  at 
10  in.  Depth  of  Penetration 
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Figure  B-53.  Lagrangian  Grid  and  Fracture  Pattern  in 

Welded  Tuff  at  15  in.  Depth  of  Penetration 


APPENDIX  C 


EXPLORATORY  PENETRATION  SOLUTION 

A short  exploratory  calculation  of  the  penetrator  firing 
at  the  Tonopah  site  was  conducted  to  assess  the  material  model, 
in  particular,  the  behavior  of  the  comminution  and  post -frac- 
ture models,  and  to  provide  tentative  pre-test  information  on 
expected  penetrator  stress  loadings  and  stress  levels  in  the 
ground.  Also,  at  this  point  in  the  program,  it  was  thought 
that  the  penetrator  firing  date  was  imminent,  and  it  became 
desirable  to  provide  some  pre-test  calculational  results.  As 
it  turned  out,  the  penetrator  firing  was  subsequently  post- 
poned, giving  time  to  determine  new  experimental  properties 
data,  revise  the  material  model,  and  perform  the  principal  pene- 
tration solution  described  in  the  main  text. 

This  solution  used  the  first  material  model  developed, 
which  included  the  newly  formulated  comminution  and  post-frac- 
ture models.  The  significant  differences  between  this  model 
and  the  subsequent  revised  model  were  that  the  comminuted  rock 
was  considerably  weaker  and  that  a lower  frictional  stress 
was  applied.  These  differences  are  illustrated  in  comparative 
plots  of  the  yield  surfaces  and  the  friction  rules.  Figures 
13  and  16,  Section  2.2. 

To  provide  a shorter  run-up  time  to  peak  loadings  on  the 
penetrator,  this  solution  was  started  with  the  penetrator  ini- 
tially buried  5 in.  into  the  ground.  The  initial  Lagrangian 
grid  is  shown  in  Figure  C-1.  The  tuff  cells  around  the  pene- 
trator nose  are  initially  unstressed  and  intact.  The  basic 
zoning  in  the  target  was  -.8  x .8  in.  cells,  or  a zone  width 
of  1/4  the  projectile  radius.  The  projectile  was  the  same  as 
previously  described  (Section  2.1),  weighing  517  lb  and  having 
an  initial  velocity  of  1500  ft/sec. 
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The  calculation  of  the  penetration  was  carried  out  to  a 
time  of  .64  msec  and  a penetration  depth  of  16.4  in.  The 
resulting  time  histories  of  penetrator  deceleration  and  axial 
forces  on  the  penetrator  are  shown  in  Figure  C-2.  These  his- 
tories have  an  interesting  character,  in  that  there  is  a rapid 
build-up  of  force  at  early  times  as  the  tuff  loads  up  to  the 
fracture  surface,  followed  by  a reduction  in  force  as  the 
fractured  material  degrades  in  strength.  After  this  initial 
excursion,  the  forces  climb  smoothly  as  the  penetrator  buries 
deeper.  This  excursion  does  not  occur  if  the  impact  starts 
on  the  ground  surface.  The  deceleration  is  seen  to  be  leveling 
off  at  the  end  of  the  solution,  having  reached  about  1800  g's. 

The  contributions  of  force  due  to  normal  and  frictional  applied 
stresses  are  separately  shown  on  the  plot.  At  the  end  of  the 
solution,  17%  of  the  total  axial  force  being  applied  was  due 
to  friction. 

The  ground  response  at  the  end  of  the  solution,  at  t=.64 
msec,  is  depicted  in  Figures  C-3  to  C-5,  which  show  the  Lagran- 
gian  grid  and  fracture  regions,  particle  velocity  field,  and 
principal  stress  field.  As  discussed  in  Section  1.3,  these 
results  indicated  that  the  material  model  appeared  to  be 
providing  a qualitatively  correct  media  response  to  penetration. 

The  normal  and  tangential  stress  distributions  applied 
along  the  penetrator  at  t=.64  msec  are  shown  in  Figure  C-6. 

The  peak  normal  stress  is  3.5  kb,  occurring  1.4  in.  back  from 
the  nose  tip.  The  peak  compressive  stress  occurring  in  the 
rock  next  to  the  penetrator  as  a function  of  depth  is  shown 
in  Figure  C-7.  Except  for  the  material  initially  next  to  the 
buried  portion  of  the  nose,  the  peak  ground  stress  has  an 
approximately  constant  level  of  3.6  kb.  Finally,  profiles  of 
radial  stress  vs  radius  along  three  horizontal  planes  next  to 
the  penetrator  are  shown  in  Figure  C-8.  These  moderate-to-shallow 
stress  gradients  appear  to  be  reasonably  well-resolved  with  the 
zoning  employed,  as  indicated  by  the  dots  on  the  curves. 
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Figure  C-3.  I.agranglan  Grid  and  Fracture  Pattern  In  Welded  Tuff,  Exploratory 
Solution 
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